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ABSTRACT 

Context. Models of many astrophysical gamma-ray sources assume they contain a homogeneous distribution of electrons that are 
injected as a power-law in energy and evolve by interacting with radiation fields, magnetic fields and particles in the source and by 
escaping. This problem is particularly complicated if the radiation fields have higher energy density than the magnetic field and are 
sufficiently energetic that inverse Compton scattering is not limited to the Thomson regime. 

Aims. We present a simple, time-dependent, semi-analytical solution of the electron kinetic equation that treats both continuous and 
impulsive injection, cooling via synchrotron and inverse Compton radiation, (taking into account Klein-Nishina effects) and energy 
dependent particle escape. We use this solution to calculate the temporal evolution of the multi-wavelength spectrum of systems where 
energetic electrons cool in intense photon fields. 

Methods. The kinetic equation for an arbitrary, time-dependent source function is solved by the method of Laplace transformations. 
Using an approximate expression for the energy loss rate that takes into account synchrotron and inverse Compton losses including 
Klein-Nishina effects for scattering off an isotropic photon field with either a power-law or black-body distribution, we find explicit 
expressions for the cooling time and escape probability of individual electrons. This enables the full, time-dependent solution to 
be reduced to a single quadrature. From the electron distribution, we then construct the time-dependent, multi-wavelength emission 
spectrum. 

Results. We compare our solutions with several limiting cases and discuss the general appearance and temporal behaviour of spectral 
features (i.e., cooling breaks, bumps etc.). As a specific example, we model the broad-band energy spectrum of the open stellar 
association Westerlund-2 at different times of its evolution, and compare it with observations. 

Conclusions. The technique we present enables simple, computationally efficient, time-dependent models of homogeneous sources 
to be constructed and compared with multi-wavelength observations. 

Key words. Gamma rays: theory - Radiation mechanisms: non-thermal - Stars: Wolf-Rayet - (ISM) cosmic rays - open clusters and 
associations: individual: Westerlund 2 



1. Introduction 

The discovery of very high energy (VHE, E y > 100 GeV) 
gamma-rays from many different astrophysical sources (for a re- 
cent review see Volk 2006) has motivated the construction of 
simple, physical models of the emission regions. An important 
ingredient in such models is the interaction of accelerated lep- 
tons with an intense ambient radiation field, especially by inverse 
Compton scattering. For example, the recent discovery of VHE 
emission from sources possibly associated with stellar clusters 
including Westerlund-2 (Wd-2) (Aharonian et al. 2007), Berk- 
87 (Abdo et al. 2007), and Cyg OB2 (Aharonian et al. 2002, 
2005a) provides evidence for particle acceleration in an envi- 
ronment dominated by hot, massive stars driving fast winds, as 
well as supernova remnants evolving in a dilute, hot medium 
(Tang & Wang 2005). The detection of VHE gamma-rays from 
high mass X-ray binary systems like LS 1+61 303 (Albert et al. 
2006), LS 5039 (Aharonian et al. 2005b), and PSR B 1259-63 
(Aharonian et al. 2005c) adds another class of galactic sources 
where inverse Compton scattering is expected to be an important 
energy loss mechanism. 

In addition to the VHE emission, these sources show ev- 
idence for non-thermal X-ray emission, which, in the case of 
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Wd-2 and Cyg OB2, is also spatially extended (Townsley et al. 
2005; Horns et al. 2007). Leptonic scenarios have been used 
both to predict and to model the multi-wavelength spectra (Kirk 
et al. 1999; Murata et al. 2003; Khangulyan et al. 2007; Hinton 
& Aharonian 2007), although hadronic models in which the 
gamma-rays are produced in nucleonic collisions also provide a 
natural explanation (Neronov & Chernyakova 2006; Horns et al. 
2007; Bednarek 2007). 

Improved modelling of these sources is clearly necessary in 
order to reach a better understanding of the mechanisms at work. 
However, even with drastic simplifications, such as the assump- 
tion of spatial homogeneity, and the ad hoc prescription of a dis- 
tribution function of injected particles, such models can be quite 
complicated to construct. Furthermore, the need to investigate a 
large volume of parameter space makes it important to find com- 
putationally efficient algorithms. 

The temporal evolution of a non-thermal electron distribu- 
tion interacting with matter and photon fields has been stud- 
ied in great detail in the past (see e.g. Ginzburg & Syrovatskii 
(1964); Felten & Morrison (1966); Blumenthal & Gould (1970); 
Mastichiadis & Kirk (1997); Sturner et al. (1997)). The problem 
can be described by an integro-differential equation including all 
relevant energy loss mechanisms for electrons (e.g. Blumenthal 
& Gould (1970); Zdziarski (1989); Coppi (1992); Kusunose & 
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Takahara (2005)) which is however numerically quite difficult to 
treat. The problem can be simplified considerably by approxi- 
mating the discrete jumps in energy suffered by an electron un- 
dergoing Compton scatterings as a continuous energy loss pro- 
cess. In this case, the kinetic equation reduces to a partial dif- 
ferential equation. This approximation, which is accurate in the 
Thomson limit, appears to be reasonably good also in the Klein- 
Nishina regime for a wide range of target photon distributions 
(see, e.g., Zdziarski 1989). 

Analytical solutions to the electron kinetic equation have 
been found for special cases e.g., when continuous or impul- 
sive injection is considered, with the electrons cooling radia- 
tively through synchrotron and inverse Compton scattering in the 
Thomson limit as well as escaping from the system (Kardashev 
1962). Although Klein-Nishina (KN) effects are well-known 
to have an interesting impact on the steady state spectra (e.g., 
Blumenthal 1971; Kirk et al. 1999; Dermer & Atoyan 2002; 
Moderski et al. 2005b) their inclusion in time-dependent mod- 
els has so far necessitated a rather elaborate numerical treatment 
(Mastichiadis & Kirk 1997; Krawczynski et al. 2002). However, 
recently, a useful approximation for the treatment of the inverse 
Compton energy losses in an isotropic photon field, including the 
transition from Thomson to the KN regime has been presented 
(Moderski et al. 2005a). In this paper, we use this approxima- 
tion to develop a semi-analytical solution for the electron kinetic 
equation. The numerical treatment of the time-dependent prob- 
lem is thus reduced to a single quadrature. We use this method 
to discuss the temporal evolution of characteristic features in the 
electron distribution as a consequence of energy losses in the KN 
regime. These effects are of interest in the sources mentioned 
above, as well as in active galactic nuclei. The paper is struc- 
tured in the following way: in §2, we derive our semi-analytical 
solution to the continuity equation describing the temporal evo- 
lution of an electron distribution suffering synchrotron, inverse 
Compton (including KN effects), and escape losses for an ar- 
bitrary, time-dependent injection spectrum. The necessary an- 
alytical calculations are summarized in Appendices A and B. 
Numerical results found using this method are presented in §3 
and a specific application to the case of the stellar association 
Wd-2 is given in §4. We close the paper in §5 with a discussion 
of the results and comment on further refinements of the calcu- 
lations. 



2. Electron energy distribution 

2.1. Formal solution of the kinetic equation 

Assume a homogeneous source into which an unspecified accel- 
eration process injects electrons at the rate Q(y, t), which sub- 
sequently cool and escape. The kinetic equation governing the 
evolution of the number N(y, t) of electrons in the source with 
Lorentz factors y in the interval dy at time t is: 

dN(y, t) d 

— jr~ = ^ {y)N{ ^ t)] ~ v esc(yW(r, + Q(y, . (D 

In writing Eq. (1) we have assumed that each individual electron 
in the source loses energy continuously, according to dy/df = —y 
(note that y is defined to be positive) and has an energy depen- 
dent probability v esc (y)df of escaping the source in a time inter- 
val df . The continuous-loss approximation is appropriate for syn- 
chrotron radiation and for Compton scattering in the Thomson 
limit, as well as for collisional losses ("ionisation losses") in 
a fully ionized plasma. In the case of Compton scattering in 



the Klein-Nishina regime, this approximation has been exam- 
ined by several authors (e.g., Blumenthal 1971; Zdziarski 1989). 
Provided one is interested only in spectral features that appear 
in the transition region between the Thomson limit and the ex- 
treme Klein-Nishina limit, it appears to be justified. In addition, 
we assume that Eq. (1) is linear in N(y, t), i.e., that the the emit- 
ted radiation has no further effect on the electron distribution, or 
on the electron injection function. This excludes situations en- 
countered in compact sources in which feedback effects can be 
important (e.g., Katarzyriski et al. 2006; Lightman & Zdziarski 
1987; Stawarz & Kirk 2007), but is not a significant restriction 
for objects such as stellar clusters. 

In the case discussed here, where the total energy-loss rate, 
y(y), and the escape rate, v esc (y), are independent of time, Eq. (1) 
can be solved by several standard methods, e.g., by means of 
the Green's function for the homogeneous equation (Ginzburg & 
Syrovatskii 1964), or using Laplace transforms (Melrose 1980). 
Following Melrose (1980), we define the quantity r(y', y) that is 
the time required for an electron to cool from a Lorentz factor of 
y' to one of y (< y'): 



T(y',y) : 



-f 

Jy 



dy" 

y(y") 



(2) 



and the quantity A(y',y) such that 1 -exp [-A (y',y)] is the prob- 
ability that an electron escapes whilst cooling from y' to y: 



rr' v 
A(y',y):= dy" ^ 

Jy J 



■Ay") 



y(y") 

This enables the solution to be written: 



N(y 



1 f°° 



-My',y) . 



(3) 



(4) 



{Q [Y,t- T(y', y)] H[t- r(y',y)] + N(y', 0)6 (t - T(y',y))} , 

where H and 5 stand for the Heaviside and delta functions, re- 
spectively. On multiplying Eq. (4) by y(y), it is straightforward 
to assign a physical interpretation to each term in this solution. 
Thus, 

y(y)N(y, t) is the number of particles per unit time that cool 
across the Lorentz factor y, 

Q-W<y)Q[y>,t -T(y',y)]H[t -r(y',y)]dy' is the rate at 
which particles injected with Lorentz factor between y' and 
y' + dy' arrive at time t at a Lorentz factor y, and 

e -<i(y ^Niy' ,Q)6(t - T(y',y))dy' is the rate at which particles 
whose initial Lorentz factor lies between y' and y' + dy' ar- 
rive at Lorentz factor y at time t . 

As we shall see below, the cooling rate is such that in the 
limit of large Lorentz factor, synchrotron losses always domi- 
nate: y ~ y 1 as y — > oo. In this case, a finite time suffices for 
an electron to cool down to y from an initially arbitrarily large 
Lorentz factor. Defining the function y\ (t) as the Lorentz factor 
at time t of a particle whose initial energy was infinite, i.e., 



t[°°,7i(0] = t 

one can rewrite Eq (4) as 



(5) 



N(y,t) 



ro 



-My',y) 



Q[y', t - T(y', y)]H[t - r(y', y)] + N(y Q , O)^™ ) 
for y <y\ (t) 



(6) 
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and 

CO 

7 

for y>y x {f) (7) 

where 70(7, t) is denned (for y < yi(f)) by the relation 

T(yo,y) = t (8) 

The exact functional form of Eqs. (6) and (7) is determined 
by the energy dependence of the radiation processes involved 
and the escape rate, as well as the injection function specified. 
As an illustration, we present the well-known case of constant 
injection of a power-law spectrum: Q(y, t) = Qoy~ p , starting at 
t = with N(y,0) = 0, where radiative cooling is due only 
to synchrotron emission: y = b s y 2 and no escape term is present 
(v esc = 0). In this case r(y', y) = {b s y'T x -{b s yT x and A(y', y) = 
0. Eq. (6) then gives (see also Melrose (1980)): 



N(y,i) = 



Qq7' p 
(P - Vb s y 1 i . 



[1 - (1 - b s yty~ l ] : Q<b s yt<\ 



b s yt > 1 



(9) 



As implied by Eq. (9) the electron spectrum steepens from y~ p 
for b s yt <sc 1 to y-^+V for b s yt as 1 with a break at y as l/b„t = 
71(f) that moves to lower energies as time progresses: the well- 
known "spectral ageing" effect. 

For an injection function which is non-zero only within a 
range of the Lorentz factor y e [y m j n , y ma x] 



Q 



Q(y, t) for y min < y 

^ 7max 



(10) 







otherwise 



it is computationally more efficient to take account of this ex- 
plicitly by modifying the integration limits in Eqs. (6) and (7). 
Specialising to N(y, 0) = 0, one finds: 

N{y,t) = jdy'e- A ^^Q[y',t-r(y',y)]H[t-T(y',y)l 



for y 2 < y < Tii 



(11) 



with N(y, t) = for y < yi and y > y max . Here, the lower and 
upper integration limits are: 

( (Train, 70) : 72 < 7 < 7min 

(r r m ,„<y f r r m ,»<r 1} = (% Jq) . ^ <y < yi (12 ) 

((7,7max) : 71 <7<ymax, 

when y m i n < yi, while 

f(ymin,yo): y2<y<yi 
(yj ,<7mm ,yl 1<ymm ) = \ (ymin,y m ax) : yi < y < y mi „ (13) 

I (y, 7max) : 7min < J < J max? 

for yi < y min . 

The value of 70 fulfill still Eq. (8), but 71 and 72 are now defined 
by: 



T(7max,7l) ~t = 0, 
T(7min,72) ~ t = 0. 



(14) 

(15) 



Eqs (14) and (15) describe the two breaks, y\ and 72 introduced 
in the evolved spectrum due to the existence in the injection 



spectrum of a minimum and a maximum Lorentz factor, 7 mm and 
7max, respectively. 72 < 71 with both of them moving to lower 
energies as the time evolves. 

As a further illustration, Appendix B presents the solution 
of Eq. (1) for a power- law type injection function with finite 
energy spectrum, where synchrotron and inverse Compton losses 
in the Thomson regime are present and electrons escape at a rate 
v e .sc 7 (see § 2.2). 

2.2. Energy loss and escape rates 

In the continuous energy-loss approximation, the total energy 
loss rate of a relativistic electron (7 » 1), averaged over an 
isotropic distribution and taking into account synchrotron cool- 
ing, inverse Compton scattering of the electrons on ambient pho- 
tons, Coulomb losses and Bremsstrahlung emission (in a fully 
ionised hydrogen gas), is given by the formula: 

7tot = b s y 2 + b ic y 2 FKx(y) 



+Mlny + b° c ) + b B y{\ny + b%). 



(16) 



The coefficients, b s , b[c, be, and bs, and the constants b° c and 
b° B are given by the following relations (Ginzburg & Syrovatskii 
1964): 



b, = — — — -wg = 1.292 x l(T 15 (fi/mG) 2 sec" 1 , 
3m e c 



(17a) 



b ic = b,— = 5.204 x l(T 2< WeV C nT 3 ) sec" 1 , (17b) 



b c = — — - = 1.491 x 10 _14 n c sec" 1 , 



4e 6 « e 
mlc 4 h 



1.37xl0- 16 n e sec" 1 , 



and 

n / mlc A \ 3 
^ =ln fe +4=^73.4, 



b% = ln2- r =0.36 



(17c) 



(17d) 



(18a) 



(18b) 



Here, cr T is the Thomson cross-section, uq, u b are the total en- 
ergy densities of the photon and magnetic fields, respectively, 
while c, m e , and e are the velocity of light in vacuum, the electron 
mass and charge, respectively. Finally, n e is the electron density 
in units of cm 4 . In the case of ionization or Bremsstrahlung 
losses (in neutral hydrogen), the corresponding terms in Eq. (16) 
should be replaced by fe c (3 In y+18.8) andfo B [ln(191)+ 1/18] * 
5.3 b B . 

The escape of particles from the source may be modeled 
in different ways. The simplest method is to assume that the 
"source" is a region of relatively strong magnetic field and high 
target photon density. When particles are carried out of this re- 
gion by the motion of the background plasma they may not only 
leave the target fields behind, but may also undergo a sudden 
drop in energy due to adiabatic expansion. This scenario implies 
that the escape rate is independent of energy, being the recip- 
rocal of the average time taken by a fluid element to cross the 
region at a given speed. 

An alternative method, (which, however, does not take ac- 
count of a possible energy loss) is to consider the particle trans- 
port as a diffusive process with the spatial diffusion coefficient 



4 



Manolakou, Horns & Kirk: Non-thermal electrons in radiation fields 



D = (Ar 2 )/2f, being either constant or a function of the energy of 
the particle or/and time, D(y, t). A common practice is to calcu- 
late the spatial diffusion coefficient in the Bohm diffusion limit, 
which is defined by setting the step size of the random walk 
of a particle to be equal to the Larmor radius, r L = ym e c/eB 
(for relativistic electrons). In this case, D(y) = r L c/3 and the 
escape rate from a spherical region with radial extension R, 
Vesc(y) - 2D(y)/R 2 is, therefore, directly proportional to the par- 
ticle's energy: 

Vesc(r) = ^scT, (19) 

where b esc = 2m e c 2 /3eBR 2 = 3.6 x 1(T 21 s" 1 (B/mG)" 1 • 
(R/pcT 2 - 

2.3. Inverse Compton scattering 

The function F KN (y) in Eq. (16) takes into account the full 
Klein-Nishina cross section for Compton scattering and is given 
by the relation: 

CO 

^kn(t) = — f /kn(7. eoK de , (20) 



where u ea is the differential photon energy density: uo - J u ea Aeo 
and eo = hv/mc 2 is the dimensionless energy of a target photon 
of frequency v. The kernel, /kn(t, eoX was first given by Jones 
(1965). It is important to note that this kernel, although originally 
derived for an isotropic field of target photons, is also valid for an 
anisotropic photon field, provided the average over an isotropic 
electron distribution is taken (see Kirk et al. 1999, Appendix A). 
A useful approximate expression in the limit y » eo, has re- 
cently been given by Moderski et al. (2005a): 

/kn(t^o) * (l+4 r e r 3/2 , (21) 
for 

yen £ 10 4 

and these authors also estimate the implications of Eq. (21) for 
the functional form of Fkn(t), f° r different photon energy dis- 
tributions. 

For a Planckian distribution of photon energies, the Compton 
energy losses for electrons are dominated by scatterings on pho- 
tons with dimensionless energy e e ff = 2.8kT/m e c 2 . For kT <sz 
ym e c 2 , the function Frn(t) can therefore be approximated by 
treating the Planckian as a mono-energetic photon distribution 
(i.e., u €o = aod(eo - e e ff))- Using this approximation and Eq. (21) 
one finds: 

Fkn(t) * (l+4yW (22) 
with 

e e ff = 2.8kT/m e c 2 

and 

B = -3/2. 

For a power-law radiation field (u 6o = aoe ff °), ^kn(t) can 
be written in the same general form as in Eq. (22), provided the 
values of the parameters e e ff and B are suitably chosen. For chq < 
-0.5 and ao > 1, inverse Compton losses in the KN limit are 
dominated by scatterings on photons with the highest (« eo,max) 
and the lowest (« eo.min) photon energies available, respectively. 
One may therefore write, e e ff = eo jmax for ao < -0.5 and e e ff = 
eo.min for ao > 1, respectively, with B = -3/2 in each case. 



E e [eV] 




,10 



T 

Fig. 1. The cooling time t coo i = yjy as a function of Lorentz factor for 
the various energy loss mechanisms for electrons including Coulomb 
scattering (dash-dotted), Bremsstrahlung (triple dotted), synchrotron 
(long dashed), inverse Compton (Thomson limit: dotted, approxima- 
tion (see Eqn.22): short dashed, no approximation: thin line), and total 
(heavy line). The parameters are n s = 1 cnT 3 , uo = 500 eV cm' 3 , 
T = 30 000 K, and B = 10 fiG, as expected for a stellar environment 
with hot young stars. 

E e [eV] 
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Fig. 2. The solution N(y, t = 10 8 yrs) using the approximate cooling 
term of Eq. (24) in an environment identical to the one used in Fig. 1 
(solid line) compared with the steady state solution using the full ex- 
pression, for F K n and including Coulomb losses (dashed line). 



In the intermediate case of -0.5 < ao < 1 the relevant pho- 
ton energy range for scattering is broader (« 1 1 Ay) and one may 
use instead the so-called "Thomson edge" or "KN cut-off" ap- 
proximation. In this range Fkn(t) is st iU given by Eq. (22) with 
e e fr = eo.max and B = a - 1 , provided eb,min/eb,max « 1- Table 1 
summarizes the values of the parameters e e ff and B for the case 
of a black-body photon targets, as well as for power-law target 
photon fields. 

The cooling timescale for electrons T coo i(y) = yjy for 
the various energy loss mechanisms, assuming an environment 
where inverse Compton cooling in the Klein-Nishina limit is 
of importance, is shown in Fig. 1. The parameters used are 
« e = 1 cm" 3 , uo = 500 eV cnT 3 , T = 30000 K, and B = 10 /zG, 
as expected in an environment dominated by hot and young stars 
which provide a relatively high temperature target photon field. 
The cooling time for Coulomb scattering and Bremsstrahlung 
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t [yrs] t [yrs] 

Fig. 3. The quantities 71 and 72 defined in Eqs. (14) and (15) as functions of time for magnetic fields and target photon fields of the same energy 
density but for two different photon temperatures: T = 30 000 K (left panel) and T = 30 K (right panel). In both panels, y s (transition from 
synchrotron to KN cooling) and 7kn (transition from KN to Thomson) are indicated by horizontal arrows. For comparison, the thin lines indicate 
the evolution of 71 and y 2 in the Thomson limit (for 72, the cooling proceeds at all times in the Thomson regime, since 7 m i„ is chosen to be 100; 
therefore, the two curves are superposed). 



losses is calculated assuming a fully ionised gas. This is a self- 
consistent assumption taking into account the high temperature 
and energy density of the ambient photon field, which is suffi- 
cient to ionize neutral hydrogen. The cooling time is dominated 
by different mechanisms in different energy ranges. At the low- 
est energies, Coulomb scattering ("ionisation loss") is the fastest 
energy loss mechanism, whereas above approximately 100 MeV, 
inverse Compton scattering in the Thomson regime takes over. 
For electron energies above a few GeV, the inverse Compton 
scattering becomes less efficient, because of the drop in the cross 
section. At energies larger than 1 TeV, inverse Compton scatter- 
ing is less efficient than synchrotron radiation, which remains 
the dominant mechanism at higher energies. Although the tran- 
sition energies described depend on the choice of the parame- 
ters (mo, « e > T, and B), the qualitative picture remains the same. 
It is interesting to compare the approximation for F kn(t) with 
the "exact" calculation shown in Fig. 1 (using the expression 
for /kn(t, fo) given by Moderski et al. (2005a)). The approx- 
imation slightly overpredicts r coo i, but is fairly accurate up to 
very high energies (~ 10% at E 10 14 eV, corresponding to 
ye e ff =! 10 4 ). However, as one can readily deduce from Fig. 1, 
the approximation for Fkn(t) becomes inaccurate only in the 
regime where synchrotron losses are likely to dominate. Even 
in binary systems, where the radiation field density can be sub- 
stantially higher than the value assumed in Fig. 1, the ratio of 
uq/ub is unlikely to be larger than (9(1000) and, therefore, elec- 
trons with ye e ff ~ 10 4 cool predominantly through synchrotron 
radiation. 



2.4. Time dependent solutions for the particle distribution 

For astrophysical environments with dilute plasma and fields of 
high energy density, the most relevant energy loss mechanisms 
for high energy electrons are synchrotron and inverse Compton 
radiation. Therefore, one may safely discard in this case the last 
two terms in Eq. (16) which then becomes: 



Photon distribution 


P 




Planck, temperature T 


-3/2 


2MT/m c c 2 


power-law, otq < -0.5 


-3/2 




power-law, -0.5 < ao < 1 


a - 1 




power-law, ao > 1 


-3/2 


P3,min 



Table 1. The values of the parameters [} and e c ff (see Eq. 22) for two 
different photon energy distributions: a (possibly diluted) Planckian of 
temperature T, or, equivalently a mono-energetic distribution with en- 
ergy e = 2.8kT '/m^c 2 , and a power-law photon field with spectral index 
cr = -dlnif^/dlneo- For the case where -0.5 < a < 1, it is fur- 
ther assumed that eomm/eomax <k 1. All approximations are valid for 
4 7 e eff < 10 4 . 



where b ls = bic/b s = mo/«b denotes the ratio of the photon field 
total energy density to the magnetic field energy density. Using 
Eqs. (22) and (23) we may write this in the general form: 



y = b s y 2 [l+b ls (l+4ye ef[ f], 



(24) 



where B and e e jf are chosen from Table 1 for the appropriate 
target photon field. Taking into account Eqs. (19) and (24), the 
functions describing cooling time and escape probability given 
in Eq. (3) can be expressed as 



(T ' r) X b s y" 2 [\+bA\+4y"e^f] 

(r ', r) = f 

J y y"[l+b is (l+4y"e eW f] 



(25) 
(26) 



y = b s y 2 [l + b is F KN (y)]. 



(23) 



where the notation b es = b &sc /b s has been used. Note that 
for y' — > 00, t(j' ,y) remains finite, whereas A{y',y) diverges. 
Physically, this is because the synchrotron cooling rate, which 
goes as y 2 dominates at high Lorentz factor, so that the time 
taken to cool to y from an arbitrarily large y' is finite. However, 
if, as in the case of Bohm diffusion, the mean free path of the 
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electrons increases in proportion to y, the probability of escape 
during this finite cooling time tends to unity as y' — > oo. 

In Appendix A, we give the analytical form of the functions 
t(j' ,y) and A(y',y) assuming that v esc (y) x y (see Eq 19). For 
a power-law type photon energy spectrum with spectral index 
-0.5 < co < 1, these functions are given in closed form only 
for the special cases ao = and ao = 1/2, although a numeri- 
cal evaluation is straightforward using Eqs (A.l) and (A. 2). The 
analytical expressions have the advantage that the solution of 
Eq. (1) is reduced to the single quadrature given in Eqs (6) and 
(7) or (11). 

In the Thomson limit for Compton scattering Fkn — L 
Therefore, assuming v esc = b esc y, one has ft = in Eqs. (25) 
and (26), and the functions r(y', y) and A(y', y) can be written 
down explicitly: 



y = (b s + b ic )y 2 = by 2 

and 



(27) 
(28) 



The solution of Eq. (1) follows from Eqs. (6) and (7). For y < 
y 1 (t) = l/(bt): 



N(y,t) 



1 
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where yo = y/(l - bty) and, for y > l/(bt): 



N(y,t) 



CO 

h^ib- fry' 



1 



1 
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(30) 



(31) 



The solution for a continuous power-law injection with finite en- 
ergy spectrum is given in Appendix B. 

In the case where the escape rate does not depend on the 
particle energy, A(y',y) = r(y' ,y)/t esc , where f esc = l/v esc = 
const. Then, in the Thomson limit of Compton scattering r(y'y) 
is given, as before, by Eq. (28), and from Eq. (6) one finds: 

7o 

N(y,t) = J-e^'^'^^IJdy'e 1 ^ 7 '^' 

7 

Q[Y, t - {byT 1 + (by'T l ]H[t - (byT 1 + (by'T 1 ]} 



+ -^e '™ ol/ ' csc jV(y ,0)e 
by 



(32) 



for y > 1 /(bt), and for y < 1 /(bt): 

N(y,t) = -}—e-' cool/ ' csc 
by 2 



J iy'e l/byk 'Q[y', t - {byT 1 + (byT 1 ] , (33) 



where the cooling time f CO oi(y) := t(oo, y) = 1 /by is the 
time taken for an electron to cool from infinite Lorentz factor 
to Lorentz factor y. For N(y, 0) = 0, with a constant mono- 
energetic injection spectrum, Q(y, t) = Q[)H( f)5{y - y,), and for 
Ji > y > yil (1 + btyi), Eq. (32) yields a steady-state solution: 



'cool /'esc + l/7i'es 



(34) 



This solution displays a cut-off towards lower Lorentz factors at 
fcooi = f esc as has been pointed out by Kardashev (1962). 

3. Numerical results 

3.1. Comparison with the steady state solution 

Ignoring escape losses, Eq. 1 in the steady state limit, dN/dt — > 
0, becomes: 



— [yN(y)) + Q(y) = 
ay 



(29) which is easily integrated to give: 



N(y) 



00 

= J Qiy'W- 



(35) 



(36) 



In this case it is straightforward to incorporate the general term 
for y given by Eq. (16). In Fig. 2, we compare this solution with 
the time dependent solution obtained by using Eqs. (6) and (7) 
with the approximate cooling term given in Eq. (24). We take 
t = 10 8 years, which is sufficiently large to allow a steady state 
to be achieved, and assume no particles are present initially, 
N(y,Q) = 0. In this figure, escape is neglected, b es = 0, and 
the ratio of the energy density of the photon field to that of the 
magnetic field is b m = 200. The remaining parameters are the 
same as in Fig. 1 Note that here, and in the following, we restrict 
our calculations to the case of a thermal distribution of target 
photons. The qualitative picture is similar to the case of a photon 
field following a power-law. 

For y ^ 100, Coulomb losses harden the spectrum to p - 1 = 
1. Since the time-dependent solution is based upon a "simpli- 
fied" energy-loss term given in Eq. (23) neglecting the Coulomb 
energy losses, the two solutions deviate in this range. Between 
y 300 andy * y KN = (4 eo) -1 ~ 18 000(7730000 K)" 1 , 
electrons cool predominantly via inverse Compton scattering in 
the Thomson limit. For values of y > yKN, the drop in the 
inverse Compton scattering cross section in the Klein-Nishina 
regime leads to a spectral hardening until synchrotron losses 
become dominant for y ^ y s = (b 2 ^ 3 - l)/(4eo) ~ 6 x 
10 5 (fc is /200) 2/3 (7730000 K)" 1 . The positions of y K N and y s are 
indicated in Fig. 2. The slight deviation observed for 7kn ^ 
y % y s is the result of the approximation used for Fkn given 
in Eqn. (21). From Fig. 2, it is obvious that (a) the approxima- 
tion used to take the energy dependence of the inverse Compton 
scattering cross section into account follows quite well the exact 
calculation and (b) neglecting Coulomb losses mainly affects the 
low energy part of the spectrum. 

3.2. Temporal evolution 

The temporal evolution of N(y, t) for an arbitrary injection func- 
tion Q(y, t) can now be calculated. In the following we inves- 
tigate the impact of the relevant parameters (b s , bi s , bes, T). In 
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order to demonstrate the effect of KN cooling on the tempo- 
ral evolution, we plot in Fig. 3 the evolution of y\ and 72 as a 
function of time for two different temperatures (T = 30000 K 
and T — 30 K) which are representative for an environment 
dominated by hot stars or by re-processed emission from warm 
dust. For comparison, the evolution of 71 and 72 in the Thomson 
limit is depicted by thin lines. In both cases, the parameters 
have been chosen similar to the ones used above (B — 10 pG, 
uq = 500 eV cirT 3 corresponding to fe; s = 200, fe es = 0) with 
Tmin = 100 and y max = 10 9 . The evolution of 72 is dominated by 
slow cooling in the Thomson limit and therefore, the two lines 
are identical. For the high temperature environment (Fig. 3a), 71 
cools via synchrotron radiation until roughly 1 Myr, when the 
energy losses in the KN limit take over. Within a few hundred 
thousand years, the transition to Thomson-limit cooling takes 
place when 71 < 7kn- Qualitatively, the evolution of 71 for a 
cooler photon field proceeds in a similar way (Fig. 3b). It is, 
however, important to note that the respective transitions occur 
much earlier. Whereas for T — 30 000 K, synchrotron cooling 
prevails for about 1 Myr, for T = 30 K, Compton cooling in the 
KN regime is reached within a few hundred years. 

The effect on the temporal evolution of a particle distribu- 
tion which is partially cooling via inverse Compton radiation in 
the KN limit, assuming continuous power-law electron injection 
with p = 2 between y m \ n = 10 2 and y max = 10 9 , neglecting es- 
cape losses (fe es = 0) and with the environmental parameters 
described above is shown in Fig. 4. The solutions N(y, t) for 
t = 10 5 , 8 X 10 5 , and 2 X 10 6 years, weighted with y 2 so that the 
injected spectrum appears constant, are compared with solutions 
found using expressions for the cooling valid in the Thomson 
limit (see Appendix B), which are depicted as thin lines. 

The solutions for t = 10 5 yrs (see also Fig. 3a) correspond 
to the case 71 > 7, (the position of 71 is marked by an arrow 
in Fig. 4). The solution taking KN effects into account shows 
a break such that N oc y - (P +1 ) for 7 > 71 while for smaller 
7kn < 7 < 7s, a slight deviation from the p — 2 injection power- 
law can be observed which marks the onset of rapid energy loss 
in the KN regime. In the solution obtained in the Thomson limit, 
the break occurs at a smaller value of 71 (see also Fig. 3a) which 
marks the transition from the uncooled part of the spectrum to 
the part which has already suffered from cooling. At the solu- 
tion at a later time (f = 8 x 10 5 yrs), when 7kn < 7i < 7s, the 
flattened spectral shape presented earlier in the steady state so- 
lution (see Fig. 2), has already fully developed as a consequence 
of the large 7 in the KN regime. The solution obtained in the 
Thomson limit and that obtained taking KN effects into account, 
agree with each other for 7 < 7kn, which marks the transition 
to cooling in the Thomson limit. Finally, for t - 2 x 10 6 yrs, 
7i < 7kn, and the further temporal evolution is well described 
in the Thomson limit. 

In the next example, we consider a higher value of 7 m ; n = 10 4 
while keeping the same values for all other parameters as in the 
examples above. Fig. 5, shows N(y, t) for three different times: 
t = 10, 10 4 , and 10 7 years. The shape of the spectrum evolves, as 
expected, mainly due to synchrotron cooling for t Sj 10 6 years. 
At later times (e.g., t = 10 7 yrs), the spectrum shows a smooth 
transition between p — 2 and p+ \ - 3 aiy ~ y s . For smaller val- 
ues of 7, the spectrum remains almost flat, showing a slight peak 
at the position of 7 m ; n . The overall shape of this evolved spec- 
trum could, in principle, mimic the case of an uncooled spectrum 
up to TeV energies (depending on the value of bQ with a soften- 
ing at 7 » 7 S resembling a much younger source. Remarkably, 
for p = 2, the effect of cooling in KN is almost hidden when 
7min ~ 7kn- For higher values of fej s , a hard particle spectrum 
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Fig. 4. The solution N(y, t) obtained including the effects of KN cool- 
ing (thick lines) and in the Thomson limit (see Appendix C) for three 
different times chosen such that 71 > y s (t = 10 s yrs), y K N < 7i < 7s 
(t = 8x 10 5 yrs), and 71 < y K N (f = 2x 10 6 yrs). The parameters used are 
identical to the ones used e.g. in Fig. 3a. The positions of y\ at different 
times as well as of 7kn and y s are indicated by arrows and vertical lines 
respectively. 

can be retained up to even higher values of 7 than those shown 
here (see also next example). 

Changing the ratio of the photon and magnetic field energy 
densities (fe; s ) at fixed photon temperature influences the loca- 
tion of the transition between inverse Compton and synchrotron 
cooling. In Fig. 6, the effect of varying fe; s from values of 100 up 
to 10 4 is demonstrated. Increasing fe; s leads to more rapid cooling 
and a more pronounced spectral hardening as a result of inverse 
Compton cooling in the KN limit. 

Finally, we consider the effect of particle escape via diffu- 
sive particle transport leaving a sphere with radius R. In Fig. 7, 
the solution N(y,t) is shown for values of fe es = 0.1, 100, and 
1000 at t — 10 6 yrs (solid lines). For comparison, the solution of 
Appendix B in the Thomson limit is also shown (dashed lines). 
All other parameters are the same as in the previous example. 
As expected, the effect of escape losses leads to a reduction of 
total particle number. For fe es = 0.1 escape losses are negligible, 
whereas for fe es = 100 and fe es = 1000, escape losses modify 
the solution considerably in comparison to the case of small fe es - 
The solution taking KN effects into account shows a stronger 
relative suppression of particle numbers for 7 > 7s than in the 
Thomson limit. This is a consequence of the fact that, in the KN 
case, radiative cooling is dominated by synchrotron losses for 
7 > 7s which implies that escape losses are relatively more im- 
portant (A.(Y, y) ~ b esc /b s ) than in the Thomson limit, where 
My' ,y) ~ besc/b with b = b s + b\c- In general, when escape 
losses are important, the spectral features characteristic for cool- 
ing in the KN limit are suppressed, and the solution resembles 
that found in the Thomson limit. Asymptotically, for very large 
fe es » fe, radiative cooling is negligible and the two solutions 
converge. 

4. Application to a stellar association: Westerlund 2 

The discovery of VHE gamma-ray emission from young stellar 
associations like Westerlund 2 (Wd-2), Cyg OB2, and Berkeley 
87 is providing first clear evidence that particle acceleration can 
occur in these systems. While the origin of the emission and the 
nature of the accelerated particles is not clear, it has been sug- 
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Fig. 5. Temporal evolution of an electron distribution injected with p = 
2 between y ra j„ = 10 4 and y raax = 10 9 in a hot (T = 30000 K) photon 
field. The evolution is shown for three distinct times: t = 10, 10 4 , and 
10 7 years. The energy density of the photon field is 200 times that of the 
magnetic field (Z? is = 200) and escape losses are neglected (fc cs = 0). 
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Fig. 6. For a similar set of parameters as described in Fig. 5 and a 
fixed time (f = 10 6 yrs) the effect of varying b is = 100, 1000, 10 4 is 
demonstrated. 



gested that gamma-rays from stellar associations could be the 
decay product of neutral mesons in inelastic scattering of nuclei 
with the ambient medium (Torres et al. 2004; Bednarek 2007) 
or due to the excitation of giant resonances of nuclei by Doppler 
boosted UV photons (Anchordoqui et al. 2007b). Here, we cal- 
culate the broad band spectral energy distribution of Wd-2 at 
different stages of its evolution. The distance to this system is 
not very well constrained (Churchwell et al. 2004). A recent es- 
timate places Wd-2 at 2.8 kpc, with an age of 2.0 ± 0.3 Myrs, 
and a total mass of 10 4 M Q (Ascenso et al. 2007). While this 
distance estimate is based mainly on near infra-red photome- 
try and colors, the spectroscopic investigation of the cluster and 
the light-curve of the binary system WR 20a has provided a 
larger distance estimate of 8.0 +1.4 kpc (Rauw et al. 2007). 
Another recent distance estimate is based upon the possible as- 
sociation of RCW 49 with a giant molecular cloud with a mass 
of 7.5 x 10 5 M Q at a kinematic distance of 6.0 ±1.0 kpc (Dame 
2007) For our calculations we adopt a distance of D — 2.8 kpc 
and note that the larger distance would affect the total energetics, 
requiring a higher injection rate. 
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Fig. 7. The effect of escape losses is demonstrated for three different 
values of £> cs = 0.1, 100, 1000 (thin, medium, and heavy lines). The 
time has been fixed to 1 Myr. The injection rate as well as the other 
parameters are identical to the cases discussed in the previous figures. 
The solid lines are obtained taking KN effects into account while for 
comparison the analytical solution for the Thomson limit is shown with 
dashed lines. 



The angular extension, 0.18 + 0.02° of the VHE gamma-ray 
emission region (Aharonian et al. 2007) implies a radial exten- 
sion of 8.8 pc • (D/2.8 kpc). The total luminosity of the member 
stars of O-type has been estimated to be Lboi ~ 3 x 10 40 ergs s _1 
(Rauw et al. 2007). However, it is well known that the number 
of UV photons emitted by the early cluster member stars is not 
sufficient to power the HII region RCW 49 in which Wd-2 is 
embedded. Based upon the estimates of Rauw et al. (2007), only 
=s 20 % of the UV photons required to power RCW 49 can be 
attributed to known O stars. Therefore, the assumed luminosity 
of the early type stars is certainly a lower limit and the total lu- 
minosity could in fact reach values beyond 10 41 ergs s _1 . It is 
noteworthy, that the population of later type stars shows a rather 
flat projected spatial distribution even beyond a distance of 4 pc 
(Ascenso et al. 2007) which would lead to a higher average en- 
ergy density of the photon field. In addition to the stellar photons, 
re-processed emission from the dust present in Wd-2 and its sur- 
rounding will provide a contribution to the seed photon field for 
inverse Compton scattering which is however not included here 
and left for a more detailed modelling. 

Taking into account just the total power emitted by the early 
type stars in the cluster, the average energy density of the hot 
photon field at a distance R corresponding to the extension of 
the VHE gamma-ray source can be estimated as: 



13 500 ■ 



eV 



Lboi 



cm 3 \ 3 x 10 40 ergs s~ 



-2 



(37) 



Based upon the extension of the VHE source of 8.8 pc, mo = 
174 eV cm 4 . However, given the presence of additional (dif- 
fuse) UV photons implied by the ionisation state of RCW 49, 
we choose uq = 500 eV cm 4 as the average energy density 
and a temperature of T - 30000 K to characterise the target 
photon field. The magnetic field is set to 10 /^G, placing it in 
approximate equipartition with the thermal energy density of 
the gas, assuming a source radius of 8.8 pc. The injection rate 
Q(y) is assumed to be constant in time with a total power of 
7 x 10 35 ergs s _1 following a power-law with p — 2 between 
ymin = 10 4 and y max = 3 x 10 8 . The required injection power 
can be compared, for example, with the total power injected by 
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Fig. 8. Broadband spectral energy distribution of Westerlund 2 (Wd 2): 
the radio point represents the integrated 843 MHz emission of the core 
given by Whiteoak & Uchida (1997) and is considered to be an up- 
per limit of a possible non-thermal component. The X-ray lower limit 
is based upon the flux observed in the spectral hard tail detected with 
Chandra (Townsley et al. 2005). We consider this detection a lower limit 
because of the limited field of view of the Chandra telescope. The upper 
limit in the X-ray band is based upon Einstein IPC observations possi- 
bly overestimating the total flux because of limited spatial resolution 
(Goldwurm et al. 1987). The VHE data is taken from Aharonian et al. 
(2007) and the GLAST sensitivity is indicated for a one year observa- 
tion. The model curves are derived for four different times assuming a 
continuous injection of a power-law with p = 2 between y min = 10 4 and 
7max = 3x 10 s in a spherical region with 8.8 pc diameter. The total power 
injected is 7 x 10 35 ergs s _I assuming a distance of D = 2.8 kpc. The 
escape loss, synchrotron, and inverse Compton cooling are taken into 
account assuming B = 10 yuG, ;/ = 500 eV cm 4 , and T = 30000 K. 
For the inverse Compton emission, the heavy line indicates the sum of 
the emissivity from the stellar light (medium curve) allowing for KN 
effects as well as from the cosmic microwave background (thin curve). 



the stellar winds driven by the early type stars (mainly the Wolf- 
Rayet stars WR 20a&b) present in Wd-2 which has been esti- 
mated to be 5 x 10 37 ergs s _1 (Rauw et al. 2007). If the emitting 
electrons are accelerated by these stellar winds, they must take 
roughly 1 % of the kinetic power of the bulk flow. 

In hadronic scenarios, on the other hand, the required effi- 
ciency is generally much larger, because of the relatively long 
cooling times involved. For example, inelastic proton-proton 
scattering producing 7r°, operates on a timescale of t(pp — > 
7T° + . . .) x 1.5 x lO 8 /^ 1 yrs (« = n\ cirT 3 denotes the ambi- 
ent medium density), which is larger than the age of the accel- 
erator, assuming this is limited to the maximum time that mas- 
sive stars drive strong and fast stellar winds in their Wolf-Rayet 
phase (this lasts approximately, twu ~ 5 x 10 5 yrs, see also be- 
low). This implies that under most conditions, the efficiency re- 
quired in hadronic models is close or larger than unity. A similar 
conclusion can be drawn for the case of gamma-ray production 
through photo-excitation of heavy nuclei, where the efficiency 
of acceleration of exclusively iron nuclei has to reach values in 
excess of 8 % (Anchordoqui et al. 2007a). 

The electron distribution N(y, t) is computed for four differ- 
ent times t = 10 3 , 10 4 , 10 5 , and 10 6 years. The synchrotron and 
inverse Compton emissivities are calculated following the ap- 
proach described in Appendix C and the resulting broad band 
spectral energy distribution (SED) is presented in Fig. 8. In 
addition to the UV photons from the early type stars, CMB 
photons are, of course, also present with an energy density 



«cmb = 0.26 eV cirT 3 . Although the energy loss due to inverse 
Compton scattering on these photons is negligible, they are nev- 
ertheless quite important for the total hard gamma-ray emission. 
This is because the inverse Compton scattering off UV photons 
is strongly suppressed, leading for sufficiently large y to emis- 
sivities which are smaller than the contribution from scattering 
with CMB photons. We therefore include in the inverse Compton 
emissivity calculation the CMB photon field and show in Fig. 8, 
the total (heavy line width) inverse Compton emission as well as 
the contribution from UV photons (medium line width) and from 
CMB photons (thin line). It is interesting to note how the inverse 
Compton emission from the CMB photons dominates for young 
sources. With increasing age, the UV-related component at about 
1 TeV increases steadily until, after * 10 6 years — roughly the 
age of this stellar cluster — it exceeds the CMB-related con- 
tribution. This effect is due to the accumulation of electrons of 
Lorentz factor (y w 10 6 ) within the source (see Fig. 5). 

The model SED is compared in Fig. 8 with VHE measure- 
ments (Aharonian et al. 2007) and with X-ray measurements 
Townsley et al. (2005) taken with Chandra (indicated as a lower 
limit) and earlier measurements taken with the Einstein X-ray 
telescope (Goldwurm et al. 1987) (indicated as an upper limit). 
Whereas the limited field of view of the Chandra observa- 
tions may underestimate the total extended X-ray emission, the 
Einstein observations may overestimate the diffuse flux, because 
of the limited spatial resolution which does not allow the sub- 
traction of the point sources present. In this figure we also show 
the total integrated radio flux of 210Jy at 843 MHz. This we 
consider as an upper limit which should not be exceeded by 
the non-thermal radio flux produced by the energetic electrons 
(Whiteoak & Uchida 1997). Overall, the model SED is in good 
agreement with the observations. The VHE energy spectrum is 
well reproduced for t > 10 5 years, and the constraints from X- 
ray and radio measurements are not violated. It is interesting to 
note that for stars with M > 25 M Q , the Wolf-Rayet phase is 
expected to last up to roughly twR ~ 5 x 10 5 years (Maeder & 
Meynet 1987), which is consistent with the injection time re- 
quired to match the data in this model. In contrast, an age of 
t > 10 5 yrs would be too long for acceleration in a supernova 
remnant, a hypothesis which also has to contend with the fact 
that, so far, no indication for the presence of a supernova rem- 
nant associated with Wd-2 has been reported. 

The predicted X-ray flux is quite close to the upper bound 
imposed by the Einstein measurements. However, it should be 
noted that this prediction is sensitive to the energy density in the 
target photon field — an increase by a factor of 2 in the energy 
density in the seed photon field (either from the late type stars 
and/or emission from dust) reduces the required injection power 
and, consequently, the predicted X-ray flux is also reduced by a 
factor of 2. We also indicate in Fig. 8 the one year flux sensitivity 
of the upcoming GLAST gamma-ray mission 1 . Our calculations 
predict that GLAST should easily detect the source. The com- 
bined energy spectrum of GLAST and H.E.S.S. would provide 
a broad band energy spectrum which in turn might allow one to 
infer, for example, the age of the source in the model scenario 
suggested here. 



1 taken from 
homepage 



www-glast.slac.stan.edu, official GLAST 
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5. Conclusions 

We have presented a method for computing the time-dependent 
electron distribution in the continuous energy loss limit for an 
arbitrary source function, under the assumptions of 

- an isotropic electron distribution 

- negligible ionization and Bremsstrahlung losses, 

- time-independent energy losses by synchrotron radiation in 
a homogeneous magnetic field and 

- by inverse Compton scattering off a spatially homogeneous 
mono-energetic or power-law energy distribution of target 
photons, 

- an escape probability proportional to the particle energy, as 
expected for diffusion in the Bohm limit. 

The key innovation of the model is that it utilises a simplified 
expression introduced by Moderski et al. (2005b) for the inverse 
Compton scattering rate (including KN-effects) to construct ex- 
plicit analytic expressions (see Appendix A) for the integrated 
cooling time {t(j' , y), Eq. 2)) and the integrated escape function 
(A(y',y), Eq. (3)). This enables the electron distribution to be 
computed in a single quadrature. 

The extension of this approach to include different en- 
ergy dependencies of the escape probability would be relatively 
straightforward. However, the treatment of multiple and/or time- 
dependent target photon fields would probably require numerical 
evaluation of the functions t and A. 

We have used this method to present an in-depth discussion 
of the properties of cooling electron distributions, and have iden- 
tified features characteristic of cooling by inverse Compton scat- 
tering: 

1. The effect of "spectral ageing" observed in the Thomson 
limit is strongly modified. While in the Thomson limit, a 
well-defined spectral break y\ changes position with time 
(cooling break), this spectral break is not evident anymore 
for y\ < y s , when KN effects are taken into account (see 
Fig. 4). 

2. For y K N < J < Js, a spectral hardening is observed as a 
result of cooling in the KN regime (see Figs. 2,4). 

3. When choosing y m ; n « y KN , the evolved spectrum shows a 
smooth transition at y s resembling a cooling break. However, 
this transition is stationary in time and may be mistakenly 
interpreted as a cooling break (see e.g. Figs. 5-6). 

4. The effect of particle escape modifies the observed energy 
spectra stronger for y > y s than for smaller values of y which 
can lead to energy spectra resembling more the ones obtained 
in the Thomson limit (see Fig. 7). 

We have illustrated our technique by applying it to high en- 
ergy observations of the open stellar cluster Westerlund 2, and 
calculating the broad band emission spectrum. Comparing our 
results with the available observations, we have identified a con- 
sistent parameters set: b is = 200, B = 10/iG, T = 30 000 K, total 
injected power in electrons 7 x 10 35 ergs s _1 , distance=2.8 kpc, 
R = 8.8 pc, y min = 10 4 , y max = 3 x 10 8 , p = 2. Although 
the uncertainty on the distance affects the values of some pa- 
rameters (eg. R, b ls , and total power), it is interesting to note 
that the suggested model of an electron accelerator embedded 
in the intense stellar radiation field which has been injecting a 
power-law type spectrum for the past few 10 5 yrs, can account 
naturally for the observed gamma-ray spectrum, without violat- 
ing the constraints derived from X-ray and radio observations. 
Furthermore, the rather long injection time favours acceleration 



in stellar winds of, for example, Wolf-Rayet stars over acceler- 
ation in supernova remnants. A leptonic scenario of gamma-ray 
production requires only a moderate fraction (~ 1 %) of the ki- 
netic power of the stellar winds to be channelled into the acceler- 
ation of a non-thermal particle population. Hadronic models, on 
the other hand, need a rather high ambient gas density in order 
to keep the required efficiency below unity. The model presented 
here can be improved further by including a more realistic mix- 
ture of different temperature photon fields with different spatial 
extensions as observed in the vicinity of Wd-2 (including a dust 
component and later- type stars present in the system). However, 
it already makes testable order of magnitude predictions of the 
flux in the range accessible to the GLAST experiment. 
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Appendix A: Characteristic cooling and escape functions 

The integrals in Eqs (25) and (26) can be rewritten as 
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^173^ (co 2 - 1) (w-V + b is ) 



(A.2) 



= 2fo es [A^g ( Vl +4e e ff7') - Ays ( y/l + 4e eff y)] 



where the indefinite integrals T^cj) and Ap(aj) can be expressed in terms of elementary functions for (negative) integer and half- 
integer values of j3. 

In the case of black-body target photons, as well as for power-law photons with index ao < -0.5, and ao > 1, one has /3 = -3/2, 
and finds, for b ls + 1 : 



r_ 3/2 (w) = a^ln^-^/V+^+a^arctan 



2o>i - b 



1/3 ^ 



1/3 



+ a^Hoj! - 1) + a® 2 (wi - lr 1 + <P vl {ph + I) -1 . 



where 



+ a®, ln(wi + ^ S /J ) + fl w /2 ln(wi + 1) 



.(4) 



(A.3) 



a ( _3/2( fe is) 



-bf{2b 2 % + fcf/ 3 - 62>f + 2^ 3 + l)/60£ _ 1)2, 



a ( _ 2 >a> is ) = -i? /3 (2# - fof/ 3 - 2£ 2/3 + 1)/ V3(6£ - l) 2 , 



-3/2 
,(3) 



',2/3/"7»,2 '1 1,8/3 



fl^'/ofe) = ' W + £° /J + 3£ 4/3 + 2b"' + l)/3(# - l) z , 

-3/2 v 1S/ is v is is is is is y ' 

a^fe) = -3b is IAQ> is - l) 2 , 
and, for Z? is = 1 : 



■2/3 . 



a ( _ 5 3 /2 (feis) = 3£ is /4fe + l) 2 , 

a ( _ 6 3 /2 fe ) = " 1/4fe + 1) ' 
flS /2 (*i.) = -l/4(*i.-l) 



r. 3/2 M = -i ln( W 2 - u + 1) + JL ln( W + 1) + 1 Info, - 1) + ^ J^^ ■ 
For the indefinite integral A_ 3 / 2 (w) one finds, for b is ± I: 

(2<o-b}P 



A_ 3/2 (w) = ln(w 2 -^ + fcf) + c P ' /2 arctan 



+ c® /2 ln(w- 1), 



1/3 



+ c®^ ln(6> + ^/ 3 ) + c% ln(w + 1) 



(A.4) 



(A.5) 



where 



c« /2 fe) = [(2^ 3 -l)^ 3 -^ 3 )]/6(^ 2 -l), 

ci 2 3 >is) = (^ /3 -^ /3 )/V3(^-l), 



'2/3 



(,4/3 _,_ ,2/3^ 



ci 4 3 ) /2 (fo is ) = -l/2(fo is -l), 
C^ /2 (*is) = 1/2(^ + 1) 



and, for b ls = 1 : 

1 , 1 

A_3/ 2 (c<j) = - ln(w -io+ 1) H arctan 

6 3 V3 



(^p) + ln ( w + 1) + \ ln ( w - 1) + ^ + • 



(A.6) 



In the case of a power-law distribution of target photons with -0.5 < ao < 1, closed form expressions of the integrals can be 
found for the special cases ao = and ao = 1/2, corresponding to (3 = - 1 and ft = -1/2, respectively. For f3 = - 1 : 



A-i(w) 



2(l+b is ) 2 



\n(oj + b ls ) + 



2(1 + bis) 



ln(o/ + bis) ■ 



2{\+hs) 2 
1 



ln(w 2 - 1) • 



1 



2(b is + 1) 



2(1 + b is ) 



ln(w z -l) . 



(A.7) 
(A.8) 
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whereas for f3 = - 1 /2 one has for + 1 : 

r_ 1/2 (w) = a {y _\ j2 ln(w + b is ) + a™ /2 ln(w + 1) + a^ /2 ln(w - 1) + a ( *l /2 (a> + l)" 1 + - I) -1 

where 

a% 2 (b is ) = ^/(fc 2 - l) 2 , a% 2 (b is ) = -l/4(fc is - 1), 



(A.9) 



a% 2 (b is ) = -bJ4(b is - l) 2 , a^ /2 (b is ) = -l/4(fc is + 1) 
a% 2 (b is ) = fo is /4(fo is + l) 2 , 

and, for b ls = 1 : 

1 1 1 w 2 -3w-2 

r_i /2 (w) = — — ln(w+ 1)+ — ln(w- 1) + —r- . 

16 16 8 (oj + l) 2 (a> - 1) 

Similarly, for the function A_i/2 we find, for b\ s + 1: 

i) 

-1/2 



A_i/ 2 (w) = c^ /2 ln(w + b is ) + c {% \ j2 ln(w + 1) + c^ /2 ln(w - 1) , 



(A. 10) 



(All) 



where 



c% 2 (bis) = bt s /0f s -l), 
c% 2 (b is ) = -l/2(b is -l) 

and, finally, for bi S = 1 : 

3 11 
A_i /2 (w) = - ln(w + 1) + - ln(w - 1) + -(co + iy l 



c% 2 (b is )=l/2(b is + l), 



(A. 12) 



Appendix B: 

For the case where accelerated particles start to inject from the central source at time r with a power-law type finite energy spectrum: 

Qoy~ p for y min < y < y max 

otherwise 



(B.l) 



undergoing synchrotron and inverse Compton losses in the Thomson regime: y = (b s + b{c)y 2 = by 2 , and escaping from the system 
(with radius R) at a rate v esc = b esc y (see § 2.2), the solution of Eq. (1) splits into two branches, depending on whether the maximum 
injected energy, y max has had enough time to "cool down" to a value y\ being smaller or larger than the minimum injected energy, 

If Tmin < 7i, then the solution reads: 

0, y < 72 

Ky-^ l \\ -[1 -by(t-toW), 



N{y,i) 



0, 

If Tmin > J\ , one gets: 
(0, 



Tmin < y < 71 
71 < 7 < 7 m ax 
7 > Tmax 



(B.2) 



#(y,0 = 



7 < 72 

K r~ s [y'rL - [y~ x - b C - Ml , 72 < 7 < n 

^T'O^L - Trntx), 71 < 7 < 7n 

Ky-^h-i^)"], 

L /max J 

0, 7 > 7n 



(B.3) 



< Y < 



TS7n 



where q = p - 1 + b esc /b, K = Qo/bq, s = 2- b esc /b and the values of 71 and 72 are given by Eqs (14) and (15), respectively. 
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Appendix C: Inverse Compton and synchrotron emissivities 

The spectrum of photons scattered by an energetic electron from an isotropic photon gas which follows a differential photon number 
density dn = n(e)de has been derived e.g. by Blumenthal & Gould (1970) in the head-on collision approximation and is given by 
their Eqn.2.48: 

dN %e 3cr T mc 3 n(e)de 

■t-t- = f(q) (Ci) 

dfdej y e 

with 

1 (4ey/mc 2 ) 2 

2 1+ Aqey/mc 

and 



f(q) = 2q In q + (1 + 2g)(l - q) + - 1 \ [[ _ / 2 (1 - q) (C.2) 



4e7/mc z (l - ei) 

The total inverse Compton spectrum from a distribution of electrons dN = N(y)dy follows from the integration over e and y: 
dN IC C C d N y6 

dj = J de J ^dTdt' (C4) 

The emitted synchrotron power per unit frequency interval emitted by a single electron with pitch angle a and magnetic field B 
is given e.g. by Ginzburg & Syrovatskii (1964): 

m = V3^sin« F /v\ (Q5) 



with the critical frequency: 

, 3eB sin a ^„ 

v c = y 2 — (C.6) 

47rmc 

and 

F(x) = x J dtK 5/3 (t). (C.7) 

X 

A simple approximation of F(x) has been used for the calculations presented in this paper (see e.g. Melrose (1980)): 

F(x) = 1.85x 1/3 exp(-x). (C.8) 

This approximation provides a relative accuracy better than 1 per cent in the region of the maximum at x ~ 0.29, and still reasonable 
accuracy in the broad range of 0.1 < x < 10. The luminosity per unit frequency interval from an electron distribution dN = N(y)dy 
is calculated by integration over y: 

J(v) = J dyN(y)j(v). (C.9) 
We assume an isotropic magnetic field such that sin a = y/2/3. 
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